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摘要 : 运动 载体 的 方位 角 信 息 对 导航 定位 十 分 重要 。 一 般 认 为 ， 全 球 定 位 系统 (Global 
Positioning System ，GPS) 等 卫星 定位 系统 无 法 利用 单 接收 机 进行 角度 测量 ， 因 此 也 就 无 法 直 
接 获 得 单 接收 机 载体 的 相关 角度 信息 。 全 球 定位 系统 等 本 身 虽 然 无 法 对 用 户 终端 进行 测 姿 ， 
但 却 能 够 在 测速 的 基础 上 解 算 获得 运动 中 接收 机 载体 的 实时 方位 角 。 对 这 一 原理 进行 了 研究 
和 阐述 ， 并 对 其 精度 水 平 进 行 了 分 析 。 通 过 基于 车 载 定位 的 实测 试验 ， 与 磁力 计 传感器 的 角 
度 测量 结果 进行 对 比 ， 对 方法 的 可 行 性 和 精度 水 平 进行 了 验证 。 结 果 表 明 ， 利 用 全 球 定位 系 
统 测速 的 方法 虽然 无 法 进行 接收 机 测 姿 方位 角 的 估算 ,但 对 接收 机 运动 方位 角 的 估算 是 切实 
可 行 的 ， 其 精度 与 载体 运动 速度 有 关 ; 当 接 收 机 载体 达到 一 定 的 运动 速度 时 ， 即 能 获得 较 高 
的 精度 水 平 。 由 此 估算 获得 的 方位 角 信 息 亦 能 作为 一 项 状态 量 作用 于 定位 结果 的 卡尔 受 滤波 
等 模型 中 ， 进 一 步 提高 滤波 模型 的 精度 。 
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N 角度 (方位 角 ) 信息 在 导航 定位 中 具有 重要 的 作用 。 在 获得 运动 载体 的 方位 角 后 ， 可 以 确定 车 辆 
的 行驶 方向 、 船 舶 的 航行 航向 ， 从 而 确定 其 运动 轨迹 。 在 车 载 导航 等 领域 ， 还 能 够 结合 电子 地 图 ， 利 
用 方位 角 信 息 更 好 地 进行 地 图 匹配 与 修正 。 并 且 还 可 以 在 对 定位 解 算数 据 进行 后 处 理 的 误差 滤波 模型 
中 引入 和 结合 角度 信息 ， 得 到 更 加 丰富 的 状态 量 ， 从 而 优化 滤波 模型 ， 以 获得 更 理想 的 滤波 效果 。 

一 般 认为 全 球 定位 系统 等 能 够 解 算得 到 用 户 的 三 维 坐标 位 置 、 三 维 运动 速度 、 加 速度 及 加 加 速度 
等 信息 ， 却 不 能 实时 进行 角度 测量 ， 也 就 无 法 直接 获得 接收 机 运动 载体 的 瞬时 方位 角 。 因 此 只 能 根据 
两 个 定位 点 进行 事后 计算 ， 以 得 到 两 点 间 的 方位 角 信 息 。 然 而 ， 这 种 做 法 缺乏 实时 性 ， 也 会 由 于 单 点 
定位 的 误差 较 大 而 导致 精度 较 差 。 为 此 ， 很 多 卫星 定位 终端 被 整合 与 安置 了 陀螺 仪 、 加 速度 计 、 磁 力 
计 等 传感器 ， 以 实时 获得 一 定 精 度 的 方位 角 信 息 。 

实际 上 ， 全 球 定位 系统 等 本 刁 虽然 无 法 对 用 户 终端 进行 测 姿 ， 但 却 能 够 测量 和 计算 得 到 接收 机 载 
体 的 运动 方位 角 信 息 。 这 使 得 在 很 多 情况 与 场合 下 ， 能 够 降低 终端 成 本 ， 不 依赖 额外 的 硬件 设备 即 能 
根据 全 球 定位 系统 自身 的 测量 数据 进行 角度 测量 。 本 文 针对 这 一 问题 进行 了 研究， 提出 了 在 全 球 定位 
系统 测速 的 基础 上 ,通过 简单 的 原理 与 方法 估算 得 到 接收 机 载体 的 运动 方位 角 信息 。 并 且 通 过 误差 传 
递 原理 对 估算 精度 进行 具体 分 析 与 预测 。 最 后 ， 在 车 载 实 测试 验 中 对 其 进行 了 分 析 和 验证 。 实 验 结 
表明 ， 基 于 简单 的 原理 与 方法 ， 即 能 在 全 球 定位 系统 速度 测量 的 基础 上 实时 估算 得 到 接收 机 运动 载体 
的 方位 角 信 息 ， 并 且 在 多 种 应 用 情况 下 ， 还 能 获得 较 高 的 精度 水 平 。 


1 卫星 导航 测速 原理 


全 球 定位 系统 接收 机 的 测速 方法 可 根据 原理 分 为 被 动 式 、 主 动 式 两 大 类 。 其 中 被 动 式 测速 应 用 较 
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为 广泛 ,易于 实现 全 天 候 、 全 测 程 和 高 精度 的 测速 "1 。 采 用 被 动 测速 方式 的 运动 载体 只 要 安装 全 球 
定位 系统 接收 机 ， 便 可 按照 用 户 所 要 求 的 时 间 间 隔 ， 在 动态 定位 的 同时 实时 测 得 整个 运行 轨道 上 各 点 
处 的 运行 速度 ， 而 运行 载体 本 身 不 需要 发 射 任何 信号 

卫星 和 用 户 接 收 机 之 间 的 相对 运动 导致 观测 到 的 卫星 信 号 频率 发 生变 化 。 这 种 多 普 蔓 频 移 通 常 通 
过 卫星 接收 机 中 的 载波 跟踪 环 测 得 。 已 知 卫 星 速度 ， 则 多 普 勒 频 移 可 用 来 计算 接收 机 速度 。 多 普 勒 频 
移 也 称 为 距离 变化 率 ， 可 视 为 相对 速度 矢量 在 卫星 视线 方向 上 的 投影 。 但 是 ， 由 于 接收 机 钟 差 变化 率 
( 即 频 率 偏 移 量 ) 的 影响 ， 测 量 结果 有 偏差 ， 实 际 测量 得 到 的 是 伪 距 变化 率 。 

由 载波 相位 测量 得 到 的 伪 距 增 量 与 平均 伪 距 变化 率 成 正比 ， 或 者 说 与 时 间 间 隔 内 用 户 相对 于 卫星 
的 速度 在 卫星 视线 方向 上 的 投影 成 正比 。 接 收 机 运动 载体 的 三 维 速度 (X, Y., 2,) 与 站 星 距 离 变 化 
A p, 之 间 的 关系 为 mE 

;2[05 - X005 - X.) + (G; - Y); - Y) + (Z; -Z20(Z; - Z0 lr; 


"a (1) 
+c(d7,—d7,) IE V5. 
其 中 ， 各 参量 的 含义 定义 如 下 : 
(1) UC. Y Z OG, Y Y;, Z;) 分 别 为 第 j 颗 卫 星 的 三 维 运 动 速度 和 位 置 ， (X,Y,, Z) 为 运动 载 
um 体 的 位 置 ， 由 全 球 定位 系统 信号 实测 结果 计算 而 得 ， 可 视 为 已 知 量 。 
(2) 为 用 户 至 第 j yey Y, 2;) 于 (X,Y ,2,) 计 算 而 得 ， 可 视 为 已 知 量 。 
(3) p; 为 站 星 距离 变化 率 ， 由 全 球 定位 系统 接收 机 按 下 式 测 得 


c p; - LN - (f, - f AT] Cf.) AT, (2) 
= 其 中 ，/ 为 接收 机 实际 接收 的 载波 频率 ; / 为 第 j 颗 卫 星 发 射 的 载波 频率 ，N 为 全 球 定位 系统 接收 机 
z 实测 的 积分 多 普 勒 频 移 计 数 ， “为 电磁 波 的 传播 速度 ，A7 为 全 球 定位 系统 接收 机 测速 更 新 率 ， 即 电 
= 用 户 拟定 的 测速 时 刻 a, ZMR. ESAE, KIE p, 也 可 视 为 已 知 量 。 

Y (4) d e, 为 接收 机 时 钟 偏差 变化 率 ( 钟 速 ) ， 一 般 为 1 ns/s。 当 测速 精度 要 求 达到 每 秒 亚 米 级 其 至 


E 更 高 时 ， 必 须 考 虑 其 影响 。 
(5) d c, 为 卫星 时 钟 偏差 变化 率 ， 小 于 0. 1 ns/s， 可 略 去 。 
(6) Li 、Li,, 分 别 为 电离 屋 、 对 流 层 的 时 延 变化 率 。 因 测速 间隔 通常 较 短 ( 秒 级 )， 故 有 i, =Li,,=0。 


(7) (X, Y, Z) 为 全 球 定位 系统 卫星 的 三 维 运行 速度 ， 有 
1 


X; -—(X - kX;) 
q 
y: = Y - kY' 
J =—( j = 2) , (3) 
q 
— 1 
=- (Z; -&Z;) 
q 


其 中 , 丰 、9 为 与 地 球 引力 、 卫 星 轨 道 等 有 关 的 参量 ， 均 可 计算 求 得 ; 上 标 1、2 分 别 表示 接收 机 在 测 
速 的 、b 时 刻 所 获得 的 数据 。 

由 上 述 对 (1) 式 的 分 析 可 见 ， 共 有 接收 机 的 三 维 速度 (X,Y,，2Z,)、 接 收 机 时 钟 钟 速 d v, 4 个 未 知 
XC. OWN 4 颗 可 视 卫 星 ， 按 照 (1) 式 联 立 4 个 方程 即 可 求解 。 在 单 点 定位 工作 中 ， 卫 星 接 收 机 即 按 
上 述 过 程 解 算 ， 获 得 测速 结果 。 由 此 还 可 进一步 求解 ， 得 到 加 速度 、 加 加 速度 等 参数 (5 。 


2 利用 速度 分 量 估算 运 动 载体 方位 角 


由 上 节 中 的 讨论 可 知 ， 利 用 4 颗 以 上 卫星 ， 可 以 解 算得 到 接收 机 载体 在 东 、 北 、 天 等 各 个 方向 上 
的 速度 分 量 。 实 际 上 ， 利 用 这 些 速度 分 量 即 可 估算 得 到 运动 载体 在 空间 的 方向 角 信息 。 特 别 地 ， 对 于 
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车 辆 、 船 舶 等 可 视 为 在 二 维 平面 上 运动 的 载体 ， 只 需 利 用 东 、 北 两 个 方向 上 的 速度 分 量 即 可 。 
如 图 1， 在 解 算得 到 接收 机 载体 在 正 东方 向 上 的 速度 向 

量 vy. 写 在 正 北 方向 上 的 速度 向 量 冰 后 ， 即 可 合成 得 到 一 个 A? 

速度 向 量 


Vi 


DA 成 =D 东 二 Dk (4) 

不 考虑 硬件 与 数据 输出 延迟 ， 在 观测 瞬间 ， 可 记 接 收 机 载体 

在 正 东 方向 上 的 行驶 距离 为 |vs | . dA7， 在 正 北方 向 上 的 行 
驶 距离 为 | | .dAT。 对 于 vs 与 正 北 方向 的 夹 角 a， 有 

v, 接收 机 


. (5) 
D 东 


由 (5) 式 可 计算 得 到 每 次 观测 瞬间 接收 机 载体 的 方位 角 a。 


图 1 利用 速度 分 量 确定 方位 角 
3 精度 总 分 析 Fig. 1 Determining azimuth angle by 


velocity component 


3.1. 影响 全 球 定 位 系统 测速 的 误差 源 分 析 

在 20 世纪 90 年 代 ， 全 球 定位 系统 接收 机 速度 估算 的 主要 误差 源 是 由 选择 可 用 性 ( Selective Availability, 
SA) 政 策 引起 的 卫星 时 钟 频率 拌 动 4 7,。 如 今 选择 可 用 性 政策 已 被 撤销 ， 该 项 误差 已 经 可 以 忽略 不 
计 。 而 电离 层 、 对 流 层 延 迟 和 多 径 效 应 在 短 时 间 内 的 变化 较 小 ， 加 上 现在 接收 机 的 更 新 频率 较 快 ， 所 
以 一 般 也 可 忽略 。 对 于 中 等 速度 以 下 的 接收 机 ， 速 度 测量 的 主要 误差 来 源 是 预测 的 卫星 位 置 (总 , Y, 
ZARE, Y, Z) 的 误差 !9 。 

从 (1) 式 可 以 看 出 ， 卫 星 轨 道 误差 对 站 星 距离 变化 率 的 影响 可 概算 为 

dp; = [UE - X,)dX, 十 m - Y, ) dY, * (Z, z Z,)aZ,]/r, ; (6) 
假定 站 星 距离 ,为 20000 km， 卫 星 运行 速度 为 3. 2 km/s。 当 卫星 轨道 误差 为 10 m 时 ， 对 站 星 距离 变 
化 率 的 最 大 影响 约 为 1.6mm/s; 当 卫 星 轨道 误差 为 100 m 时 ， 其 影响 增加 到 1. 6 cm/s。 目 前 ， 国 际 全 
球 卫 星 导航 系统 服务 ( International Global Navigation Satellite System Service, IGS) 公 布 的 全 球 定位 系统 
广播 星 历 精度 为 1.6 m， 了 卫星 钟 差 的 精度 为 7 ns'”"。 因 此 ， 卫 星 轨 道 误差 对 全 球 定位 系统 测速 精度 的 
影响 很 小 。 

对 于 卫星 速度 误差 而 言 ， 则 通过 方向 余弦 直接 作用 于 站 星 距离 变化 率 。 在 取消 选择 可 用 性 政策 
后 ， 利 用 广播 星 历 计 算 卫 星 速度 精度 优 于 1 mm/s 。 

如 果 用 户 是 高 动态 的 ， 那么 还 有 男 外 的 问题 。 由 于 距离 增 量 只 给 出 一 个 时 间 间 隔 内 的 平均 速度 ， 
所 以 很 显然 ， 高 加 速度 和 速度 突变 的 情况 会 带 来 误差 。 另 一 方面 ， 这 种 高 动态 还 会 使 得 接收 机 相位 品 
声 变 大 ， 从 而 引入 额外 的 误差 。 由 于 这 些 误差 的 存在 ， 美 国 国防 部 并 没有 给 出 全 球 定位 系统 速度 估算 
的 具体 性 能 指标 "“; 。 根 据 经 验 来 看 ， 现 在 市 场 上 大 多 数 的 全 球 定位 系统 单 频 导航 接收 机 ， 其 三 维 方 
向 上 的 测速 精度 一 般 在 0.3 m/s 以 内 (10o); 而 基于 差分 全 球 定 位 系统 (Differential Global Positioning 
System, DGPS) 、 载 波 相 位 差分 技术 (Real-Time Kinematic, RTK) 等 方法 的 测速 精度 ， 则 能 够 达到 厘米 
/ 秒 甚至 更 高 5 。 

3.2 方位 角 估 算 精 度 分 析 

由 (5) 式 可 知 ， 在 二 维 平面 上 获得 方位 角 的 精度 主要 取决 于 接收 机 在 正 东 与 正 北 两 个 方向 上 的 速 

EDE vg, vo Wlogl=x, |vg |=y， 对 (5) 式 进行 偏 微分 处 理 可 得 


da = (arctan z) |xdx + (arctan 2) | ydy 


F! 
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由 误差 协 方差 传播 率 !2 可 得 方位 角 a 的 方差 为 
n een (8) 
X ty X ty i 
为 简化 公式 ， 假 定 接收 机 载体 在 正 东 与 正 北方 向 上 的 速度 误差 大 小 相同 ， 有 o,-7o,-c,, W8) 
式 可 转化 为 


c, c, 
au = iy me (9) 
x xy lex | + [vi | 
其 中 ，wV |v# l+ v | 即 为 载体 在 二 维 平面 上 的 合成 运动 速度 |vaw |。 故 (9) 式 可 进一步 简化 为 
E (10) 
Dv 


由 (10) 式 可 以 看 出 ,估算 得 到 的 方位 角 精 度 取决 于 载体 在 二 维 平面 上 速度 的 大 小 | vss | 以 及 载体 
在 正 北 和 正 东 方向 上 速度 的 误差 o,。 由 上 文 可 知 ， 对 于 一 般 的 全 球 定位 系统 单 频 导航 接收 机 ， 其 e， 
的 大 小 一 般 不 超过 0. 3 m/s(10o); 而 基于 差分 全 球 定位 系统 、 载 波 相位 差分 技术 等 方法 的 测速 精度 ， 
则 能 够 达到 厘米 / 秒 甚至 更 高 。 因 此 ， 不 妨 分 别 取 o,= 0.3 m/s 及 o,=0.05m/s， 并 根据 |v | 的 值 对 
c, 的 估算 精度 进行 考察 。 


R1 方位 角 估 算 精 度 预测 


Table 1 Prediction of estimated azimuth's accuracy 


[van | (ms) 1 5 10 20 50 

(。) 0,=0.3 m/s 17. 19 3.44 1.72 0. 86 0. 34 
Oo 

i c, 20.05 m/s 2. 86 0. 57 0.29 0. 14 0. 06 


从 表 1 中 可 以 看 出 ， 当 该 方法 应 用 在 行人 的 运动 方向 判定 时 ， 由 于 vag | 的 值 较 小 (人均 步 行 速 
度 约 为 1~2 m/s) ， 因 此 方位 角 佑 算 的 误差 偏 大 。 而 当 应 用 在 车 载 导 航 情 况 下 时 ， 由 于 车 辆 的 平均 行 
驶 速度 约 在 10-20 m/s 之 间 ， 所 以 能 够 获得 较 精确 的 方位 角 佑 算 值 : 对 于 一 般 单 频 单 点 定位 的 全 球 
定位 系统 民用 接收 机 ， 其 方位 角 佑 算 误 差 优 于 3?"(1z) ; 而 基于 载波 相位 差分 技术 等 高 精度 定位 方法 
的 方位 角 佑 算 精 度 ， 则 优 于 0.5"(lc)。 


4 ”实测 试验 及 结果 分 析 


为 验证 方法 的 可 行 性 与 有 效 性 ， 于 2013 年 10 月 在 北京 进行 了 车 载 定 位 实验 。 将 一 个 典型 磁力 计 
模块 固定 安置 在 车 辆 项 盖 上 ， 且 其 正 北 指向 与 车 辆 车 头 方向 平行 ， 如 图 2。 


图 2 磁力 计 安 置 示意 图 


Fig.2 A magnetometer mounted on the car 
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首先 ， 将 磁力 计 模 块 固定 不 动 ， 以 测量 其 真实 角度 测量 误差 。 数 据 输出 窗口 如 图 3， 其 中 最 右 列 
为 磁力 计 和 角度 测 量 值 ， 输 出 间隔 为 0. 5 s。 经 过 一 小 时 的 静止 测量 ， 计 算得 到 该 磁力 计 模 块 的 角度 测 
量 误差 为 1 r=0.88"。 另 外 , 在 实验 中 还 需 注 意 对 磁力 计 的 输出 数据 进行 磁 偏 角 的 修正 。 在 北京 地 
区 ， 磁 偏 角 值 约 为 5. 833 ” 。 
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图 3 


磁力 计 角 度数 据 输出 窗口 


Fig.3 The output window of magnetometer angle data 


在 完成 上 述 准 备 工 作 后 ， 车 辆 开始 分 别 以 约 5 m/s, 10 m/s, 20 m/s 的 速度 行驶 大 致 相等 的 时 间 。 
与 此 同时 ， 通 过 串口 实时 读 取 一 典型 全 球 定 位 系统 单 频 接收 机 定位 模块 的 输出 数据 (包括 单 频 接收 机 
获得 的 星 历 、 伪 距 等 信息 ) ， 利 用 上 述 方法 解 算得 到 车 辆 行驶 的 瞬时 方位 角 ， 并 与 磁力 计 和 输出 的 方位 
角 信 息 进行 比较 ， 两 者 之 间 的 差异 如 图 4。 
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图 4 磁力 计 输 出 的 方位 角 与 解 算 所 得 方位 角 的 比较 


Fig.4 Comparison of calculated azimuth angle and azimuth angle output by the magnetometer 
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需要 对 图 4 说明 的 是 ， 实 际 上 解 算得 到 的 方位 角 与 磁力 计 输 出 的 方位 角 之 间 总 是 存在 大 约 5° 的 
固定 偏差 。 这 主要 由 采用 的 坐标 系 与 磁力 计 的 坐标 系 的 偏差 引起 。 男 外 ,磁力 计 本 身 的 系统 测量 偏 
差 ， 模 块 固定 在 车 辆 上 时 与 车 辆 前 后 轴 方 向 上 的 偏差 以 及 磁 偏 角 修 正 的 残 差 也 会 引入 系统 误差 中 。 这 
个 固定 偏差 不 会 影响 实验 结果 。 图 4 中 的 统计 结果 是 由 实测 结果 减 去 固定 偏差 后 得 到 的 。 

记 由 本 方法 解 算得 到 的 方位 角 为 ga, ， 人 磁力 计 输 出 的 方位 角 为 a,， 两 者 之 间 的 差异 为 Aa: 


Aa =Q; -o,. (11) 
根据 图 4 中 的 数据 ， 由 (12) 式 计算 得 到 Aa 的 均 方 根 误差 ru : 
es. AES; (Aa)! , (12) 
其 中 , n 为 观测 历 元 的 个 数 。 而 根据 误差 传播 率 有 
gu PU nO (13) 
可 得 
gu OA. T, : (14) 


又 已 知 ww=0.88 ， 可 由 (13) 式 计算 得 到 各 不 同行 驶 速度 下 的 O, o 
按 上 述 过 程 分 别 计算 车 辆 以 5 m/s, 10 m/s, 


20 m/s 的 速度 行驶 时 ， 得 到 的 方位 角 误 差 vc。 。 ee 
Ri ci 结果 见 表 2 Table 2 Azimuth calculation error with 
ü zu (un o 
PP different velocit 
- 可 以 看 出 ,， 表 2 中 的 实验 结果 与 表 1 中 的 预 一 
um 测 分 析 基 本 一 致 ， 验 证 了 方法 的 可 行 性 和 精度 。 "en | m/s DE 
Oa 4. 45? .32? . 19? 
5 2 二 i e 5 2.32 1. 19 


由 于 全 球 定位 系统 本 身 角 度 及 姿态 测量 能 力 的 不 足 ， 为 实现 对 接收 机 终端 的 角度 测量 ， 在 很 多 定 
位 终端 被 整合 与 安置 了 陀螺 仪 、 加 速度 计 、 人 磁力 计 等 传 感 混 。 实 际 上 ， 全 球 定位 系统 等 本 身 虽 然 无 法 
对 用 户 终端 进行 测 姿 以 获得 各 种 姿态 角 ， 但 却 能 够 在 测速 的 基础 上 ， 实 时 解 算 获 得 接收 机 载体 的 运动 
方位 角 。 本 文 对 这 一 原理 进行 了 研究 和 阐述 ， 并 对 其 精度 水 平 与 误差 指标 进行 了 推导 和 分 析 。 最 后 ， 
在 车 载 定 位 的 基础 上 ， 通 过 与 磁力 计 和 角度 测量 值 的 对 比 ， 对 该 方法 进行 了 验证 。 

实验 与 分 析 结 果 表明 ， 利 用 单 频 全 球 定位 系统 测速 进行 接收 机 运动 方位 角 估算 的 方法 切实 可 行 ， 
且 原理 简单 、 计 算 方 便 。 方 法 的 精度 与 接收 机 载体 的 运动 速度 有 关 ， 当 载体 达到 一 定 的 运动 速度 时 ， 
即 能 获得 较 高 的 精度 水 平 ， 能 够 满足 一 般 性 导航 定位 的 精度 需求 。 根 据 这 个 特性 ， 该 方法 尤其 适用 于 
车 载 定位 、 船 舶 导航 等 领域 ， 能 够 不 依赖 于 磁力 计 、 陀 螺 仪 等 外 部 设备 获得 一 定 的 方位 角 测 量 精 度 ， 
具有 良好 的 应 用 价值 。 

另外 值得 指出 的 是 ， 由 于 能 够 获得 较 高 的 精度 ， 因 此 该 方法 估算 得 到 的 方位 角 信 息 还 能 作为 一 项 
状态 量 作用 于 对 定位 解 算 结果 进行 后 处 理 的 卡尔 曼 滤 波 等 模型 中 ， 从 而 优化 滤波 方法 ， 获 得 更 理想 的 
滤波 效果 。 这 也 是 值得 深入 研究 和 探讨 的 一 个 方向 。 
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Azimuth Estimation of a Moving Receiver by GPS Speed Measurement 
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Abstract: The azimuth of a moving receiver is of vital importance for navigation and positioning. 
Generally, satellite systems such as GPS are considered that they cannot measure azimuthal angle with a single 
receiver, so, related angle information of a single receiver cannot be obtained directly, either. This paper 
proposes that these systems can get the real-time azimuth of a moving receiver through real-time resolution on 
the basis of speed measurement, albeit they themselves may not be able to measure the posture of the user 
terminal. This paper also studies and elaborates the principle above and analyzes its precision. The feasibility 
and precision of this method are verified by comparing the angle measurement results of magnetometer sensor in 
the actual measurement experiment of vehicle positioning. Results show that using GPS speed measurement to 
estimate the azimuth of a moving receiver is feasible, although it is impossible to estimate the azimuth of the 
posture of the user terminal. And the precision is related to the speed of the receiver—when the receiver moves 
up to certain speed, high precision can be obtained. Moreover, the azimuth information gotten by the method 
can also be used as a quantity of state in such models as Kalman filtering of positioning results and further 
improve its accuracy. 
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